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Fracture Analysis of Notched 

(*) 

Composites 

by 

F. Delale ' ' 

Drexel University, Philadelphia, PA 19104 


Abstract 

In this study the problem of part-through cracks emanating 
from a central notch in a composite laminate is considered. The 
composite laminate is modeled as an orthotropic plate and along 
the branches, it is assumed that the stresses. are carried by the 
fibers only. The stresses along the uncracked portions are 
assumed to be proportional to the crack opening displacement at 
that point. After formulating the branched crack problem, the 
model is applied to a (±45°) 2 s composite laminate. The stress 
concentration factors for the fibers are computed for . various 
values of the spring constant K. Sample results showing the 
crack opening displacement are also displayed in the figures. 


1 . INTRODUCTION 

The load carrying capacity of a composite laminate may be 
reduced significantly due to notches or crack— like defects. 

Thus, notched composites have been studied extensively. In cer- 
tain composites, cracks making an angle with the pre-existing 
notch may emanate from its tip. These cracks are part-through, 
meaning that only some of the laminae are broken and they tend 
to meander along the fracture path [1,2]. As the geometry sug- 
gests, an exact stress analysis of the problem is extremely com- 
plicated. In this study, the following approximate model is 
used in the analysis. The laminate is assumed to be an ortho- 
tropic plate. The stresses in the unbroken laminae of the 
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part-through cracks (along the fracture path) are represented by 
springs and at each point are assumed to be proportional to the 
crack opening displacement, i.e., c = KAv. 

Starting with the formulation of the crack geometry shown in 
figure 2, the through crack problem is formulated in terms of the 
displacement derivatives of the crack surfaces . This formulation 
yields a system of singular integral equations which can be 
solved numerically. Next, the part-through crack problem is formu- 
lated, by relating the crack surface tension to the crack opening 
displacement. A new set of singular integral equations is obtained. 
At the crack tips the stress intensity factors are computed. Also 
the stress concentration factors at the tip of the original slit 
are computed for various values of K. These factors can be used 
in studying fiber damage. Also some examples showing the crack 
opening along the fracture path are displayed in the figures. 


2. FORMULATION OF THE PROBLEM 

To formulate the part- through crack problem, we start with the 
geometry shown in figure 2. In this symmetric configuration, all 
the cracks are assumed to be through cracks. Since the problem is 
linear, the solution will be constructed by the superposition of 
several cracked plates. 

2. a. The Centrally Cracked Plate 

Consider an infinite orthotropic plate with a central crack 
extending from -a to +a. The governing equation of the problem 
can easily be written in terms of the stress function F(x,y) as 
follows [3,4]: 
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E , E » G and v being the engineering elastic constants for 
x ' y ' xy xy 

the orthotropic material. 

Define: 


oo — S rv v 

F (x,y) = / g(a,y) e da 


(4) 


Substituting in (2) we obtain: 

B - B a 2 + a 4 g = 0 

1 dy q z dy z 

The solution of (5) may be written as: 
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where and co 2 are the roots of 
B^ 4 - B 2 to 2 + 1 = 0 


(7) 


such that Re(w^) > 0 and Re{u> 2 ) > 0. 

Taking the symmetry of the problem into account, one may con- 
sider one half of the plane only and equation (6) may then be 
rewritten as: 

_UJ | CX J y _ _ ^ 2 | a |y 
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Thus, the stress function F becomes: 


F(x,y) = / [Ae 


, | -tD. a y -rax 

w ll a l y + Be 2 ] e da y>0 (9) 



4 . 


Taking derivatives of the stress function, 
following expressions for the stresses: 
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To formulate the crack problem, the 
taken as the unknown function. 
Defining , 
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and using the following stress-strain relationship, 
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Integrating (13) with respect to y, and then differentiating the 
resulting expression of v(x,y) with respect to x, we have: 
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Noting that: 


iim 9v = f( . 
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( 15 ) 


equation (14) yields: 


f(x) = i /“ A | a | a [(a 12 to 1 - ^) - (a^u^-a^ ~1 } ] e ~ iaXd “ ( 16 > 
-°° 1 to_ 


The inverse transform of (16) gives: 


2 

oj 1 t0 ? 


A (a) = 


2ir^ a | a 


>2 2 . — a 

a 22 {u) l” w 2 


f f(t)e lat dt 


(17) 


Finally, substituting the expression of A (a) into expressions 
(10) and integrating over a the stresses can be written in terms 
of the unknown function f(t) as follows: 
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The expressions given in ( 18 - 20 ) may be used to obtain the stresses 
at any point in the plane in terms of the unknown function f(t). 


2 .b. Formulation of the inclined crack 

Consider now a crack making an angle 0 with the x-axis. For 
this case the problem will be formulated in the s-n coordinate 
system. Again referring to [ 3 , 4 ] the governing equation for an 
anisotropic material in terms of the stress function F 2 (s,n) can 
be written as: 
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To simplify the formulation, we note that for (± 45 °) composites 
the angle 0 defining the path of fracture is close to 45 ° and 
E x = Ey (see for example [ 1 ] ) . Thus, for 0 = 45 ° and E x = E^ we 
have: 
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and equation (21) reduces to: 
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Again we use the Fourier transform to formulate the problem. 
Defining 


F 2 (s,n) = /" h (n r a) e“ ias da 


(26) 


substituting into (24) and after some algebra we obtain: 
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— 00 


In these expressions r^ and r 2 are the roots of: 


4 2 
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such that Re(r^) > 0, Re(r 2 ) > 0. 

Using the continuity conditions along the axis n 
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Then, the expressions given in (27) become: 
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Since the inclined crack under consideration is denoted by A in 
figure 2, from here on, for the quantities related to this crack, 
a superscript A - will be used. 

Dif ferentiating (33) , the stress expressions are found to be 
for n>0 , 
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Again, the problem will be formulated in terms of the crack sur- 
face derivatives. Thus, defining: 
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Substituting (43) and (44) into expressions (35-40) and performing 
the infinite integrals over a, we obtain: 
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For cracks B/ C, D shown in figure 2 , similar expressions may be 
obtained. Taking into account the symmetry conditions, for each 
crack the stresses may be written as: 
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3. THE INTEGRAL EQUATIONS AND THE RESULTS 
3 .a. The case of through cracks 

For the through crack configuration shown in figure 2, the 
integral equations can be obtained by superimposing the solutions 
for the central crack and the inclined cracks A, B, C, D. The 
infinite plate is under a uniform tension c parallel to the y-direc 
tion. For the uncracked plate, the stresses along the- crack lines 


A 


are: 
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along y = 0 , a = a , t = o 
y o xy 

2 

along n = 0 , a n = a Q cos 0 = a Q / 2 


T ns = °o sin0cose = a 0 //2 


Thus, the crack surface tractions may be written as: 
a (x,0) = -a n 


°n (s ' 0) = ~ a o^ 2 
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With the use of equations (18-20) and (48-59), the boundary 
conditions given in (62-64) yield the following system of singular 
integral equations: 
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for n = 0 
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To complete the formulation of the problem, one has also to 
write the following single-valuedness conditions for the dis- 
placements : 


/ a f (t)dt = 0 (68) 

-a 


f C 
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g^tjdt = 0 


(69) 


/ c g„ (t)dt = 0 (70) 

b ^ 

The singular integral equations (65-67) may be solved by one of 
the collocation methods given in [5-7] . However, to apply the 
methods, first the system must be normalized. Thus, we make use 
of the following transformations: 

for equation (65) , 

x = az -a<x, t<a 

t = ap -l<z, p<l 

f(t) = f*(p) 


( 71 ) 
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for equations (66,67), 
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By substituting (71) and (72) into (65-70), the system will be 
normalized. In this case, since only through cracks are con- 
sidered, the important fracture parameters, i.e., the stress 
intensity factors are computed. The stress intensity factors are 
defined as: 
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Noting that the normalized unknown functions will have the 
usual square-root singularity at both ends, one may write: 
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(80) 


where the functions F(p), G-^t^ and G 2 (t 1 ) are bounded in the 
interval [-1,1]. Thus, with the definitions given in (73-77), we 
obtain: 


k l (c) = " ft 2b il^ r 2 +r l )r l r 2 G2(1) 

k l = </r 2b ll ( r 2 +r i ) r i r 2 


(81) 

(82) 

(83) 


k 2 (c) 
k 2 (b) 


^ 2b 1 J(r 2 +r 1 ) G l (1) 

= & 2b 11 (r 2 +r 1 ) G l (_1) 
with £ = (c-b)/2 


(84) 

(85) 


Results have been obtained for a (-45) 2s composite. The 
material properties used in the computations and given below were 
taken from [8] . 


E x = 

19.6 

GPa 

E y = 

G 

19.6 

GPa 

33.7 

GPa 

xy 


v , = 

0.735 


yx 




Some sample results are presented in tables (1,2) and figure 3. 
First, the results were checked for the case when the inclined 
cracks were located far away from the central crack. For the nor- 
malized stress intensity factor at the tip of the central crack. 
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the expected value of 1 was recovered. Also the normalized stress 
intensity factors at the tips of the inclined cracks reached the 
expected limiting value of 0.5. Table 1 shows the results for an 
inclined crack whose center is kept fixed, but the length of which 
varies. As expected the stress intensity factors at the tips 
x = +a and s = b, assume higher values as the two cracks approach 
each other. Table 2 and figure 3 show the results for a cracked 
plate, where the length of the inclined cracks is kept constant. 

Again, as expected, the stress intensity factors at tips x = +a, 
and s = b decrease as the cracks are further apart from each 
other. 

3 .b. The branched crack case 

In this case we still have through cracks, but we assume that 
all cracks are joined. Thus, we obtain a centrally cracked plate 
with branches running at ±45° from the tips of the central crack. 

The formulation of the problem is the same except that one has to 
put b = 0 for the limits of the integrals in equations (65-67) and 
the single valuedness conditions (68-70) are not necessary anymore. 
However, by putting b = 0, the nature of the singularity at x = a 
and s = 0 will change since the integral equations will now have 
generalized Cauchy kernels. A brief inspection of the crack 
geometry shows that, the crack displacement derivatives must be 
bounded at x = a, and s = 0, since all wedge angles are less than 
it. (For stress singularities in isotropic wedges see [9] and for 
stress singularities in orthotropic wedges see [10,11].) The solu- 
tion of singular integral equations with generalized Cauchy kernels 
may be obtained by the methods suggested in [12-13] . In this study 
the system is solved by assuming that the crack surface displacements 
for the central crack are symmetric and that the crack surface deriva- 
tives are bounded at x=a and s=0 . Thus, the conditions (68-70) are 
replaced by: 

f 1 f* ( p ) dp = 0 (86) 

-1 

G 1 (-l) = 0 (87) 

and 

g 2 (-d = 


0 


( 88 ) 
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The numerical computations show that, the convergence of the 
system is excellent and as expected F(l) = 0 since the crack sur- 
face displacement derivative is known to be bounded at x = a (or 
2 = 1) . For this case the only relevant fracture parameters are 
the stress intensity factors at the tip of the branched crack at 
s = +c. The normalized stress intensity factors k^(c) and ^(c) 
are given in table 3 and figure 4. The results indicate that the 
normalized stress intensity factors decrease with increasing crack 
length. 

In literature, only for the isotropic case one may find some 
results for the branched crack. Therefore, it is very difficult 
to compare the results directly. As a check, the anisotropic 
results found in this study were compared with the isotropic 
results found in [14] . Although the material and geometry were 
close but not exactly the same, the results were found to be of 
the same order. 

3.c.' A central crack with part-through branches 

The foregoing formulation was performed mainly to study this 
case which appears in (±45) and other composites. In a centrally 
cracked plate under uniaxial tension, some inclined cracks which 
make up an angle with the initial central crack may develop. How- 
ever, some of the laminae are not broken, i.e. , the cracks are not 
through cracks. In this work, the following model is adopted: 

It is assumed that the unbroken laminae are weakened and they 
carry only tension or compression. Furthermore, it is assumed that 
the unknown normal stress in the unbroken laminae is proportional 
to the crack opening displacement at that point, i.e.: 

a — KAv (89) 

where K is called the spring coefficient and Av is the opening 
of the crack. Thus at each point of the fracture path one may 
compute a bounded stress and thus a stress concentration factor. 

It is expected that generally the stress concentration factor will 
be larger near the original crack tip, namely, near s = 0. The 
spring coefficient in relation (89) may be a constant giving a 
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linear stress-crack opening relationship or may vary from point 
to point making relation (89) non-linear. This relationship 
which must be determined experimentally, in general is believed to 
be non-linear [15] . However, in this study, for simplicity it is 
assumed to be linear, thus K is taken as a constant. 

To incorporate this relationship in the formulation, we may 
replace boundary condition (63) by: 

a n (s,0) =-Oq/ 2 + KAv (90) 

or 

a n (s,0) = - a Q /2 - K / c g 2 (t)dt (91) 

s 

The problem may then be solved as in case (b) with the difference 
that, in the singular equations system (65-67) , equation (66) must 
be replaced with equation (91) . After solving the system, one may 
compute the stress intensity factors k^(c) and k 2 (c) at the tip of 
the branched crack, the normal stress distribution a along the 
fracture path, the stress concentration factor for the fibers at 
s = 0 (and x = a) . As stated earlier, the spring coefficient K is 
determined experimentally. In [15] the value of K is given for a 
certain composite. In this study, results are obtained for four 
different values of K, but chosen of the same order of K given in 
[15]. The stress intensity factors are given in table 4, and the 
stress concentration factors are presented in table 5 and figure 5. 
The results indicate that the stress concentration factors which 
may be used to study the breaking of fibers near the tip of the 
original crack, vary significantly with the spring coefficient K 
and the ratio c/a. As K increases, i.e., for stiff er materials, 
the stress concentration will be higher, meaning that for such 
materials it is more likely that the fibers ahead of the slit will 
break. Here it must be noted that, the stress concentration factor 
(SCF) is not always maximum at s = 0. In some instances, the maxi- 
mum stress concentration factor was found to be maximum at a small 
distance from the tip of the original slit. Two examples showing 
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the crack opening along the fracture path are displayed in figures 
6 and 7 . Figure 6 shows the crack opening for a short crack 
(c/a = 1) and figure 7 shows the crack opening for a larger crack 
(c/a = 5) . It is observed that the crack profiles may be signi- 
ficantly different and for the long crack the maximum stress con- 

'V 

centration occurs not at s = 0 but at s/c = 0.2. 
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Table 1. Stress Intensity factors for the cracked plate 
(through cracks) 



k x (a) 

k x (c) 

k 1 (b) 

k 2 (c) 

k 2 (b) 

b/a 


o 0 VT 

°(A 



mwm 

1.003 

0.535 

0.537 

0.577 

0.580 

115 

1.075 

0.532 

0.544 

0.570 

0.592 

1.25 

1.175 

0.531 

0.548 

0.565 

0.602 

1.0 

1.332 

0.531 

0.550 

0.560 

0.617 

0.75 

1.571 

0.532 

0.549 

0.554 

0.641 

0.50 

1.955 

0.536 

0.545 

0.548 

0.691 

0.25 

2.717 

0.545 

0.563 

0.548 

0.823 

O 

i — 1 
• 

o 

4.003 

0.560 

0.677 

0.559 

1.098 


Table 2. Stress Intensity factors for the cracked plate 
(through cracks) 

U = = 1) 



k x (a) 

k 1 (c) 

k 1 (b) 

k 2 (c) 

k 2 (b) 

b/a 

a 0 V^ 

°o^ 

° ( y r 


a of 

0.1 

2.745 

0.593 

0.728 

0.597 

1.019 

0.5 

1.607 

0.548 

0.582 

0.574 

0.689 

1.0 

1.332 

0.531 

0.550 

0.560 

0.617 

1.5 

1.221 

0.523 

0.534 

0.548 

0.584 

2.0 

1.161 

0.518 

0.526 

0.540 

0.564 

2.5 

1.125 

0.514 

0.520 

0.533 

0.550 

3.0 

1.100 

0.512 

0.516 

0.527 

0.541 

o 

• 

in 

1.050 

0.507 

■ 0.509 

0.515 

0.520 

10.0 

1.017 

0.503 

0.503 

0.505 

0.506 



































































































26 . 


Table 3. The stress intensity factors for a completely cracked 
laminate 

{Branched crack case) 





























Table 4. The stress intensity factors at the tip of the branched crack. (Part-through cracks case) 



K = 0.5 x 10 7 l/m 

K = 10 7 1/m 

K = 5 x 

10 7 1/m 

K = 10 8 1 

/m 

c/a 

k l ^ 

k 2 (c) / a 0 '/cT 

k i /O ofi 

k 2 (c)/o 0 fe- 

k l ^ /a 0 ^ 

k 2 (c) /o 0 £" 

k l /o 0^ 

k 2 (c) /o Q /c 

0.25 

1.211 

0.725 

IBBB9I 

0.741 

0.765 

0.820 

0.567 

0.863 

0.50 

0.820 

0.662 

0.729 

0.679 

0.417 

0.739 

0.295 

0.762 

1.0 

0.563 

0.610 

0.469 

0.624 

0.236 

0.660 

0.165 

0.669 

2.0 

0.394 

0.561 

0.308 

0.571 

0.144 

0.587 

0.100 

0.590 

5.0 

0.248 

0.509 

0.182 

0.512 

0.081 

0.516 

0.055 

0.516 
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Table 5. The stress concentration factors for the stressed 
fibers at the tip of the notch 



K=0.5xl0 7 1/m 

K=1Q 7 1/m 

K=5xl0 7 1/m 

K=10 8 1/m 

c/a 

SCF 

SCF 

SCF 

SCF 

0.25 

0.19 

0.35 

1.19 

1.75 

0.50 

0.24 

0.42 

1.12 

1.50 

1.0 

1.28 

0.45 

0.93 

1.12 












Figure 3. Variation of the normalized stress intensity factors 
with b/a for a cracked plate U = D • 
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Figure 5. The stress concentration factors for various values 
of the spring constant K. 


Avx 10 7 (rn) 


O.G I 




Figure 6. The crack 
K = 10 7 1/m. 


opening displacement Av for c/a 
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Figure 7. The crack opening displacement Av for c/a = 5, 
K = 10 8 1/m. 
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